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The effective linear and third-order nonlinear optical properties of metallic colloidal crystal immersed in a 
graded-index host fluid are investigated theoretically. The local electric fields are extracted self-consistently 
based on the layer-to-layer interactions, which are readily given by the Lekner summation method. The 
resultant optical absorption and nonlinearity enhancement show a series of sharp peaks, which merge in 
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^ . a broadened resonant band. The sharp peaks become a continuous band for increasing packing density 

Q , and number of layers. We believe that the sharp peaks arise from the in-plane dipolar interactions and the 

' surface plasmon resonance, whereas the continuous band is due to the presence of the gradient in the host 

refractive index. These results have not been observed in homogeneous and randomly-dispersed colloids, 
, and thus would be of great interest in optical nanomaterial engineering. 
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There have been considerable interests in searching for materials with a large nonlinear sus- 
ceptibility. The most common way to achieve this is to use composite materials in which the 
constituent components possess large intrinsic nonlinear responses. Noble metal (typically gold, 
silver and copper) is often chosen as an ingredient due to their extremely large and fast nonlinear 
optical response. Many different microstructures have been exploited in an attempt to access the 
intrinsic optical nonlinearity of metals, for example, in random metallodielectric composites 
fractal films,- and alternative bilayers,^ etc. They basically rely on the enhanced local fields in 
space or on the effectively lengthened scale of the interactions between the matter and the inci- 
dent light field. However, there is also a great demand for particular optical materials in devices 
applications, which would benefit from additional tunability of the optical properties. Recently, 
we studied graded composites, which provide an extra degree of freedom for controlling the opti- 
cal properties of these materials."^ In fact, there exist in nature abundant graded materials, such as 
biological cells^ and liquid crystal droplets.- Furthermore, many artificially-graded-index optical 
metamaterials and elements have been fabricated nowaday s.- 

CoUoidal crystal has been extensively studied in nanomaterials engineering and its potential 
applications range from nanophotonics to chemistry and biomedicine.- It is desirable and of inter- 
est to use dielectric-coated metallic nanoparticles with varying shell thickness to form a dielectric 
constant gradient.- In this Letter, we theoretically investigate a metallic colloidal crystal immersed 
in a graded-index host and demonstrate a giant enhanced optical nonlinearity band, which is con- 
trollable by the gradient and by the easily-tunable colloidal structure as well. These colloidal crys- 
tals can be prepared via layer-by-layer self-assembly, templated sedimentation, methods based 
on capillary forces, and electrokinetics. ■^'^'iSiiJJ^ We basically use the quasistatic point-dipole ap- 
proximation, which suffices in terms of characterizing both the gradient effects and the lattice 
effects, otherwise the solution is formidable, either from a Green's function formalism or from 
first-principles .— 

The theoretical calculations are deployed on a model tetragonal lattice with uniaxial 
anisotropyii (see Fig. [T]). As examples, let us focus on the body-centered-tetragonal (bet), the 
body-centered-cubic (bcc) and the face-centered-cubic (fee) lattices respectively. Among the three 
cases bcc lattice has the lowest packing density while fee lattice has the highest one.-^ Extensions 
to other colloidal structures such as simple tetragonal lattices are straightforward and similar re- 
sults are expected. Taking advantages of the normalized interlayer interaction tensor T (i.e., Tij 
denotes the interaction strength between two in-plane dipole arrays) given by the Lekner summa- 
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tion method,— we solve the following self-consistent equation for the layer dependent local field 



where a is the lattice constant (b and c are the other two lattice constants) as shown in Fig.[T], aj 
is the layer-dependent linear bare polarizability, here i, j label the crystal layer and denotes the 
total layer number. Note E.^-* in the self-consistent equation is not simply the externally applied 
electric field Eq due to the presence of gradient. However, it is the field inside the graded host 
medium, which is thus determined by virtue of the continuity of the normal component of the 
electric displacement D in the longitudinal case (i.e., Eq parallel to the uniaxial axis). It is the 
2;-axis as shown in Fig. [l] in our case. Nevertheless in the transverse case (i.e., Eq perpendicular 
to the uniaxial axis), we exactly used the applied field Eq because the boundary condition now be- 
comes the continuity of the tangential component of electric field. We shall compare the effective 
linear and nonlinear optical responses of colloidal crystals with the different lattice structures (e.g., 
bet, bcc, and fee), made of metallic nanoparticles of linear dielectric constant ei and third-order 
nonlinear susceptibility xi, suspended in a host fluid of (see Fig.[l]). Both the longitudinal (L) 
and transverse (T) results will be discussed. The self-consistent equations over i = 1,2 - ■ ■ , N are 
then combined together to take into account the lattice effect and are being able to be transformed 
into a matrix form as E=TAE+E^°-'. More precisely, in the longitudinal and the transverse cases, 
E = {E^^''^^ } is simply N-dimensional vector and A is x diagonal matrix of the polarizability, 
which relates the induced dipole moment Pj of the particle and the local field Ej in the layer ii, and 
in fact consists of isotropic linear and nonlinear contributions. That is Pj = OjEj + /?j|EipEj/3, 
where = emr^i<^i — em)/(ei + Se^). Here r is the radii of the metallic nanoparticles. In the 
case of weak nonlinearity in the colloidal particles, i.e., xi |Ej p ^ ei in the nonlinear relationship 
Dj = eiEj + xi|EjpEj, by the perturbation expansion method,'^ one has 



It is noteworthy that the linear local fields Ej around the particles in the layer £i are actually 
obtained by assuming no intrinsic nonlinear response, i.e., we set Xi = for solving the self- 
consistent equations [Eq. dU)], which is appropriate provided that the nonlinear responses are much 
less than the linear ones. Next we use the resultant linear local fields Ej to extract the enhancement 
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factor of the effective nonlinear susceptibility Xi^^ 
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where p is the filling fraction of the metallic colloidal nanoparticles. If the nanoparticles are touch- 
ing, i.e., with the geometric constraint a^+6^+c^ = 16r^ (see Fig. [T]), then p = 7r[(g'^ + 2)/g]^/^/24 
where q = (c/ a)^/^ quantifies the degree of anisotropy of the periodic lattice - and also determines 
the normalized interlayer interaction T. This results in structure-controllable optical properties. 
Note that the average (...) in Eq. © is taken over the layers ii (i = 1, ■ ■ ■ , N) instead of over 
the nanoparticles spatial volume, because in the dipole approximation the local fields inside each 
of the particles are homogeneous. We also assume no nonlinear response in the host, which is in 
fact relatively neglectable comparing to that in the metal. Additionally, a gradient of the dielectric 
constant of the host fluid is introduced along the uniaxial direction of the colloidal crystal, i.e, 
Cm = em(-2i) in our case, where Zi represents the 2;-coordinate of layer £j. In this regard, we treat 
the host as a continuously-layered film, thus explicitly has E'-°^ = {E^^^ = Eo/em{zi)} in the lon- 
gitudinal case. The formation of the gradient in the host might be achieved by dispersing different 
polymers in it, by selectively filling with microfluidic materials,'^ or induced by the presence of 
a temperature gradient, etc. One can also simply coat the nanoparticles with different coverage 
shells. But it still remains a challenge because the novel properties from our prediction require a 
reasonably large gradient in the dielectric constant of the host. 

Figure |2l shows in logarithmic scale the optical absorption, i.e., the imaginary part of the ef- 
fective linear dielectric function Im(eeft) (upper panels) and the modulus of the nonlinearity en- 
hancement factor 7 (lower panels) defined in Eq. (jSj), as functions of the reduced frequency. We 
specifically compare the results of bet (q = 0.87358), bcc (q = 1.0), and fee (q = 1.25992) as 
shown in the three columns for = 25 layers. The dielectric function of the metallic colloidal 
nanoparticles is simply adopted as the Drude form ei = 1 — ujp/{uj'^ + iujT), while the graded di- 
electric constant in the host medium is assumed as em(2;j) = 1. + 1. 252/A^ fori = 1,2,---A^. The 
presence of the inhomogeneity in the host fluid obviously leads to a broadened and giant enhanced 
resonant band (solid lines) in the low-frequency region. This is interesting for potential telecom- 
munication applications. The results of the same colloid crytal suspended in homogeneous host 
medium with = 1-0 (dotted lines) or = 2.25 (dashed lines) are also presented, in an attempt 
to demonstrate that the broadened resonant band in some sense stems from the hybridization of 
the non-graded structures. From the absorption spectrum and the enhancement in the third-order 
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nonlinear susceptibility, we would expect an attractive figure of merit."^ That is, the materials ef- 
fectively exhibit large nonlinearity and relatively small absorption. This is certainly superior to 
pure metal because it generally has large nonlinearity and unwanted absorption concomitantly. In 
fact, the optical absorption arises from the surface plasmon resonance, which is obtained from the 
imaginary part of the effective linear dielectric constant e^f! that is extracted from the generalized 
Clausius-Mossotti formula.-i^ Note that the plasmon resonant peaks in the cases of homogeneous 
host with periodic arrangment of metallic particles, i.e., the peaks in the dashed and dotted curves 
in Fig.|21are red-shifted with respect to the corresponding ones (not shown) of a homogeneous host 
with randomly-dispersed nanoparticles predicted by the Maxwell-Garnett theory. We actually set 
a gap (e.g., coated layer thickness in experiments) of 2a/5 between the nearest in-plane lattice 
particles in order to avoid severe complications arising from the multiple image interactions. 
Thus we have the metal filling fraction p = {3/5fn[{q^ + 2)/qf/^/2A,e.g., 0.15080, 0.14692, and 
0.15994 for the bet, bcc, and fee lattices, respectively. The introduction of this gap indeed makes 
the nanoparticles size, and thus the dipole factor relatively small and somewhat unfavorably sup- 
presses the effect arising from variation of the lattice structure, as seen in Fig|2l 

Furthermore, due to the fact that we treat the continuous variations of dielectric function in the 
host as layered ones in obtaining and Pi, and the fact that the dipoles are actually distributed in 
discrete lattice nodes, a series of sharp peaks are also observable in Fig.|2l The peaks are merged 
in the broadened band and become more notable for increased gradient in the host dielectric con- 
stant (not shown), whereas they tend to disappear as the number of crystal layers is increased. 
This is also understandable in the generalized Bergman-Milton spectral representation in graded 
composites.^" In detail, the merging of the peaks into a continuous band is explicitly shown in 
Fig. |3l where we increased the layer to = 50. The fee lattice is taken as an example and we 
present both the longitudinal (left panel) and the transverse (right panel) results. The peaks in 
Fig.|2tc) and Fig.|2tf) are distinctly smeared out in Fig.|3ta) and Fig.|3tb), respectively. The trans- 
verse results [see Fig. ^c) and (d)] in the presence of the gradient is slightly different to that in 
the longitudinal case [Fig.|3ta) and (b)], but still retain the broadened bands. We ascribe this to 
the fact that the layer-to-layer interactions fall off exponentially due to the screening effect in the 
lattice,^^ therefore give no much layer- structure-dependent difference in the two cases. Note that 
the longitudinal and transverse results of crystals in homogeneous host (dotted and dashed lines in 
Fig.|3l) do not differ much as well. 

In conclusion, we theoretically exploit the optical resonant enhancement due to lattice effect 
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and gradient effect in colloidal crystals, which are made out of suspended metallic nanoparticles in 
a graded-index host. The gradient in the fluid and the colloid structure are easily subjected to tun- 
ability, for example, the structure transformation might be induced by electrorheological effects or 
by self-assembly of two kinds of particles with biochemically different surface properties, etc. In 
addition, one can also use metal-covered mangetic nanoparticles and control the suspension struc- 
ture by external magnetic field, consequently realizing magneto-controlled optical properties.— In 
this case, the electro-magnetorheological effects is a good candidate as well. Devices that could 
benefit from these materials include optical switches, optical limiters, as well as biosensors, etc. 

This work was supported by the RGC Earmarked Grant. We thank Dr. J. P. Huang for critical 
reading of the manuscript. 
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Figure Captions 



FIG. 1 : Schematic diagram of the colloidal crystal immersed in a graded host medium. In the numerical 
calculations a = bis assumed, which forms a square dipole lattice in xy-plane, while c is the lattice constant 
in the direction of z-axis. 

FIG. 2: The effective linear response Im(eeff) and third-order nonlinearity enhancement factor |7[ of the 
periodic colloidal nanoparticles immersed in a graded host fluid (solid line), and in homogeneous host 
fluids with dielectric constant £2 = 2.25 (dashed line) and €2 = 1.0 (dotted line), respectively; (a), (d) 
for bet, (b), (e) for bcc; and (c), (f) for fee, respectively. Parameters: F = 0.02a)p, = 25, a = 1, and 
|Eo| = 1. 

FIG. 3: Same as Fig.|2j but with totally = 50 layers for fee lattice, (a) and (b) the longitudinal case; (c) 
and (d) the transverse case. 
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